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Abstract 

In this paper we study the asymptotic behavior of a family of poly- 
nomials which are orthogonal with respect to an exponential weight 
on certain contours of the complex plane. The zeros of these polyno- 
mials are the nodes for complex Gaussian quadrature of an oscillatory 
integral on the real axis with a high order stationary point, and their 
limit distribution is also analyzed. We show that the zeros accumulate 
along a contour in the complex plane that has the S'-property in an 
external field. In addition, the strong asymptotics of the orthogonal 
polynomials is obtained by applying the nonlinear Deift-Zhou steepest 
descent method to the corresponding Riemann-Hilbert problem. 

1 Introduction 

1.1 Oscillatory integrals 

We study the limiting behavior of the zeros of the polynomials that are 
orthogonal with respect to an oscillatory weight function of exponential type 
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along a path T in the complex plane, 



/ n n (z)z k e lzr dz = 0, 



k = 0, 1, . . . , n — 1, 



(1.1) 



./r 



such that the integral is well defined. Parameter r > 2 is an integer and we 
will focus mainly on the case r = 3. 

Our motivation originates in a Fourier-type integral on a finite interval 
of the real axis of the general form 



where u > is a frequency parameter, / is called the amplitude and g is the 
phase or oscillator. Integrals of this type appear in many scientific disciplines 
involving wave phenomena, such as acoustics, electromagnetics and optics 
(see for example jH] and references therein). For u ^> 1, integrals of this 
kind are a recurring topic in asymptotic analysis, and we recall in particular 
the classical method of steepest descent, which can be applied when / and g 
are analytic in a neighbourhood of [a, b], see for instance |27J. 

We will concentrate on the case where the oscillator g has a single station- 
ary point £ of order r — 1 inside the interval [a, b], with r > 2, i.e., g^'(^) = 0, 
j = 1, . . . , r — 1 but g( r '(£) 7^ 0. Without loss of generality, we take this point 
£ to be the origin and the canonical example is the following: 



with a < and b > 0. Assuming a single stationary point, the general 
form ( II. 2p can always be brought into this form by a change of variables. 

1.2 Numerical evaluation 

We assume / analytic in a complex neighbourhood of the interval [a, b\. As 
shown in [TJ], one possible numerical strategy for the evaluation of (II. 3p is 
to consider paths of steepest descent stemming from the endpoints and from 
the stationary point. In this way, we can decompose the original integral as 
follows: 




(1.2) 
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Figure 1: Approximate contours of integration in the complex plane corre- 
sponding to even r (left) and odd r (right). 



where the paths are depicted in Fig. [TJ Making an appropriate change of 
variables, the line integrals along these paths have the form 



with fi = 1 for the endpoints and \i = r for a stationary point of order 
r — 1 . Each of these integrals can be efficiently approximated using Gaussian 
quadrature, because the optimal polynomial order of Gaussian quadrature 
translates into optimal asymptotic order in this setting: for n quadrature 

2n + l 

points the error behaves like 0{u f ) as u — > oo, see |8j. This order is 
approximately twice that of a classical asymptotic expansion truncated after 
n terms. 

There are two paths of steepest descent originating from the stationary 
point, called Tq and Tq in Fig. [TJ Both paths are straight lines and their 
structure depends essentially on the parity of r: for odd r, these lines form 
an angle equal to n — ir/r, whereas for even r they form one straight line in 
the complex plane. 

In order to keep the total number of function evaluations to a minimum, 
it is desirable to evaluate both line integrals using only one quadrature rule 
of Gaussian type [8] . This amounts to constructing a quadrature rule for the 
functional 





(1.4) 
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Figure 2: Location of the quadrature nodes for r = 3 on [—1,1], correspond- 
ing to n = 10 (left), n = 20 (center) and n = 40 (right). In dashed line, the 
paths of steepest descent from the origin. 

where r = Tq U Iq is the concatenation of the two steepest descent paths 
through the origin. 

In the case r = 2, this leads to classical Gauss-Hermite quadrature, which 
involves the weight function e~ x on the real line (—00, 00). Higher even val- 
ues of r lead to straightforward generalizations, and in all cases the quadra- 
ture points lie on the paths of steepest descent. 

For odd r, the functional f ll.4|) is indefinite and the existence of orthogonal 
polynomials is not guaranteed a priori. Nevertheless, the orthogonal polyno- 
mials and their zeros can be computed numerically. However, one finds that 
the zeros, which are the complex quadrature points for the integral (3), do 
not lie on the paths of steepest descent anymore. Instead, they seem to lie 
on a curve in a sector of the complex plane bounded by the paths of steepest 
descent. Their location for r = 3 is shown in Fig. [2] for several values of n. 
Similar phenomena are observed for larger odd values of r. 

1.3 Orthogonality in the complex plane 

The problem of Gaussian quadrature leads to the study of polynomials n n (z) 
that are orthogonal in the sense of (II. ip . where r is a positive integer (r > 3 
for a non-classical case) and T is the combination of two paths of steepest 
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descent of the exponential function e %%r from the origin, so 

71 2mTT 

&rgz = 1 , m = 0, 1, . . . , r — 1. 

2r r 

The straight lines Tq and in Fig. [2] correspond respectively to m = and 
m = [|J- In the case r = 3, these lines form angles of n/6 and 5tt/6 with 
respect to the positive real axis. 
Putting A n = (n/r) 1 ^ and 

P n {z) = \- n ir n {\ n z) (1.5) 

we note that (11.11) can be written in the form 

J P n {z)z k e~ nV{z) dz = 0, k = 0, . . . , n - 1, (1.6) 

where 

^(z) = -z^/r. (1.7) 

Note that (jl.5p is again a monic polynomial, and the zeros of P n (z) and 7T n (z) 
are the same but for rescaling with the parameter A n . 

The orthogonality (I1.6P is an example of non-Hermitian orthogonality 
with respect to a varying weight on a curve in the complex plane. A basic 
observation is that the path T of the integral in (11.61) can be deformed into 
any other curve that is homotopic to it in the finite plane, and that connects 
the same two sectors at infinity. For any such deformed T we still have the 
orthogonality condition (11.61) . 

In order to find where the zeros of P n {z) lie for large n, we have to select 
the 'right' contour. Stahl [25] and Gonchar-Rakhmanov [131 Sec. 3] studied 
and solved this problem, and from their works it is known that the appropri- 
ate contour should have a symmetry property (the so-called S-property) in 
the sense of logarithmic potential theory with external fields. We recall this 
concept in the next subsection. 

In the case (11.71) with r = 3, we can identify the curve with the S- 
property explicitly as a critical trajectory of a quadratic differential. Other 
cases where the potential problem is explicitly solved include [13] , [21] and [2] 
in connection with best rational approximation of e~ x on [0, oo), [J] in con- 
nection with a last passage percolation problem, and [18], [19], [17], [22], [23] 
in connection with classical orthogonal polynomials (Laguerre and Jacobi) 
with non-standard parameters. 
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Varying orthogonality on complex curves is treated in detail in the more 
recent accounts [3] and [7] , which contain a Riemann-Hilbert steepest descent 
analysis in a fairly general setting, assuming the knowledge of the curve with 
the ^-property See also [6] for an approach based on algebraic geometry and 
Boutroux curves and [5] for extensions to L p optimal polynomials. 



1.4 The S property 

Let V be a polynomial. We consider a smooth curve T C C, such that 
the integral in (11.61) is well-defined, and we want to minimize the weighted 
energy: 

I v {y) = J J log 1 - dv{z)dv{s) +Re J V(s)du(s), (1.8) 

among all Borel probability measures v supported on T. Following the gen- 
eral theory of logarithmic potential theory with external fields, see [21] , this 
problem has a unique solution, which is called the equilibrium measure on 
T in the presence of the external field Re V. We denote this equilibrium 
measure by \x. 
Let 

U»{z) = [ log— ^—dfx(s) (1.9) 



;i.io) 



\z — s\ 

be the logarithmic potential of \i. It satisfies 

2U' x (z) +ReV(z) =£, ^Gsupp/i, 
2U fl {z) + ReV{z) >£, zer\supp/x, 



for some constant £, see [23] . If T is an analytic contour, then supp /x will 
consist of a finite union of analytic arcs. Now we can define the ^-property 

Definition 1.1. The analytic contour T has the S -property in the external 
field Re V if for every z in the interior of the analytic arcs that constitute 
supp fi, we have 

[2U»{z)+ReV{z)] = -— [2U^(z) +ReV{z)]. (1.11) 



dn+ dri- 

A. 

dn± 



Here -Tp— denote the two normal derivatives taken on either side ofT. 
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The result of Gonchar-Rakhmanov then reads (for the special case of 
polynomial V): 

Theorem 1.2. Gonchar-Rakhmanov Sec. 3] If T is a contour 

with the S -property (11. lip in the external field KeV , then the equilibrium 
measure fi onT in the external field Re V is the weak limit of the normalized 
zero counting measures of the polynomials P n defined by the orthogonality 

(HU). 

1.5 Outline of the paper 

In the next section we present the main results of this paper, corresponding 
to fll.7p with r = 3. These can be summarized in the following points: 

• We present a finite curve 7 C C, which is a critical trajectory of a 
certain quadratic differential Q(z)dz 2 , see Theorem 12.11 

• We prove that this curve 7 can be prolonged to 00 in a suitable way, 
thus obtaining a curve T with the S"-property in the presence of the 
external field Re V, see Theorem 12.21 

• As a consequence of Gonchar-Rakhmanov theorem, it is possible to 
obtain the weak limit distribution of the zeros of P n {z) as n — > 00, see 
Theorem 12.31 

• Additionally, a full Riemann-Hilbert analysis of this problem is feasi- 
ble and yields both existence of the sequence of orthogonal polynomials 
P n {z) f° r large enough n and the asymptotic behavior of P n {z) in var- 
ious regions of the complex plane as n — > 00, see Theorem 12.41 

2 Statement of results 
2.1 Definition of the curve 7 

In the case f ll.7j) with r = 3, the curve with the S'-property is given in terms 
of the critical trajectory of the quadratic differential Q(z)dz 2 , where 

Q(z) = --(z + i) 2 {z 2 -2iz-3). (2.1) 
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Figure 3: The contour T consists of the critical trajectory 7 and its analytic 
continuations 71 and 72. 

The polynomial (12. ip has a double root at z = — i and two simple roots 
at z\ = — v2 + i and Z2 = \/2 + i. 

The critical trajectory 7 is an analytic arc from z\ to Z2 so that 

— / Q 1/2 (s)ds (2.2) 
m J Zl 

is real for z G 7, see [26]. We first show that this curve indeed exists. 

Theorem 2.1. There exists a critical trajectory 7 of the quadratic differential 
Q(z)dz 2 , where Q(z) is given in (12. ip . £/ia£ connects the two zeros z\ = 
—\[2 + % and z 2 = y/2 + i of Q. 

The proof of the theorem is contained in Section [3l 
2.2 Contour with 5-property 

In what follows we use the analytic arc 7, whose existence is guaranteed by 
Theorem 12.11 with an orientation so that z\ is the starting point of 7 and 
Z2 is the ending point. The + side (— side) of 7 is on the left (right) as we 
traverse 7 according to its orientation, as shown in Figure [3j 

From now on the square root Q 1//2 (z) is defined with a branch cut along 
7 and so that 

Q 1/2 (z) = -^ 2 -- z + °(^) (2-3) 

as z —7- 00. This branch is then used for example in (12. 2p . We use Q + (s), 
when s G 7, to denote the limiting value of Q l l 2 (z) as z approaches s G 7 
from the + side. 
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The curve 7 has an analytic extension to an unbounded oriented contour 

T = 7l U 7 U 7 2 (2.4) 
that we use for the orthogonality (11.61) . The parts 71 and 72 are such that 

4>\{z) = / Q 1 ^ 2 (s)ds is real and positive for z G 71, (2.5) 

and z 

4>2{z) = / Q 1 ^ 2 (s)ds is real and positive for z G 72- (2.6) 



The main result of the paper is then the following. 

Theorem 2.2. The contour T is a curve with the S -property in the external 
field ReV^. In addition, we have that the equilibrium measure on T in the 
external field Re V is given by the probability measure 

d^(s) = —Q 1 ' 2 (s)ds. (2.7) 

The proof of this theorem is presented in Section HI The general result of 
Gonchar-Rakhmanov, see Theorem II .2[ then implies: 

Theorem 2.3. Assume r = 3. For large enough n the monic polynomial 
P n (z) of degree n satisfying the orthogonality relation (11.61) with V(z) given 
by (11 -7p exists uniquely. Furthermore, denoting by 

~L,n; • • • j Z n ^ n 

its n zeros in the complex plane, we have 

(a) as n —¥ 00 the zeros accumulate on 7; 

(b) the normalized zero counting measures have a weak limit 



1 n 

where \i is given by (12 



n 

k=l 
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2.3 Riemann-Hilbert analysis 



Theorem 12.31 also follows from a steepest descent analysis for the Riemann- 
Hilbert problem that characterizes the polynomials P n . We include an expo- 
sition of this method in Section [5] for two main reasons: 

• The Riemann-Hilbert analysis not only provides the limit behavior of 
the distribution of zeros of P n (z), but also strong asymptotics of the 
orthogonal polynomials in the complex plane. 

• The present case can be viewed as a simple model problem for Riemann- 
Hilbert analysis in the complex plane. We hope that the present anal- 
ysis can be useful as an introduction to this powerful method. 

The Riemann-Hilbert problem for orthogonal polynomials was found by 
Fokas, Its and Kitaev [12], and the steepest descent analysis of Riemann- 
Hilbert problems is due to Deift and Zhou [TT] . The steepest descent analysis 
for orthogonal polynomials with varying weights on the real line is due to 
Deift et al, see [TO] and [9]. 

The extension of this method to orthogonal polynomials on curves in the 
complex plane is not new. It has already been presented in various papers, 
see for example [2] , [3] , [7] , [TO] and [TS] . However, an attractive feature of the 
example treated here is that all quantities in the analysis can be computed 
explicitly. In that respect it is similar to [18J. 

In order to formulate the additional asymptotic results that follow from 
the steepest descent analysis, we introduce some more notation. We use the 
function 2 (z) defined in f)2.6p and the related function (the g- function) 



in terms of the equilibrium measure /j on 7. In the case r = 3 there is an 
explicit expression for 2 (%)'■ 



g( Z ) = -V(z)-(j) 2 {z)-l, I 



3 + 2 l0g2 ' 



(2.8) 



which has an alternative expression 




z(z + z 2 — 2iz — 3 



log(z - i + Vz 2 - 2iz - 3) + - log 2. (2.9) 



10 



The so-called global parametrix N(z) is defined in terms of the function 

with the branch cut taken along 7. N(z) is a 2 x 2 matrix valued function 
with entries 

iVn(z) = n^) = M±^r!, WuW = ^ = «*) - «*)- 



2i 



that also appear in the asymptotic formulas in Theorem 12.41 
Finally, near the endpoint Z2 we require a conformal map 



/(*) 



2/3 



(2.11) 



which maps 7 and 72 near Z2 into the real line. A local Riemann-Hilbert 
problem is solvable explicitly in terms of the usual Airy function Ai(z) and 
its derivative Ai'(z). These functions appear in the asymptotic formula in 
part (c) of Theorem 12.41 that is valid in a neighborhood of z 2 . 

The steepest descent analysis of the Riemann-Hilbert problem then leads 
to the following result: 

Theorem 2.4. Assume r = 3. LetUg(zi) andUsfa) be small neighbourhoods 
of the points z\ and Z2 given before. As n — > 00 the polynomial P n (z) has the 
following asymptotic behavior: 

(a) Uniformly for z in compact subsets of C \ 7, we have 

P n (z) = N u (z)e n ^ (1 + 0(l/n)) (2.12) 

as n —7- 00; 

(b) There is a neighbourhood U of 7 in the complex plane, so that, uni- 
formly forzeU\ (U s (z 1 ) U U s (z 2 )): 

P n (z) = e n [^-'] (e- n ^ z) N u (z) ± ie n<t>2{z) N l2 (z) + 0(l/n)) , (2.13) 

where the + (—) sign in (12.131) is valid for z in the part of U that lies 
above (below) the curve 7; 
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(c) Uniformly for z G Us(z 2 ), we have as n —¥ oo, 

p n ( z ) = y?e n [^"'] (n 1/6 f 1/ \z)/3' 1 (z)Ai(n 2 / 3 f(z)) (1 + 0{l/n)) 

- n -V*f-W( z )P( z )M\nWf(z)) (1 + 0(l/n))) , 

with the same constant I as in (12. 8p . 



3 Proof of Theorem 2.1 



It follows from the general theory, see that three trajectories of the 
quadratic differential Q(z)dz 2 emanate from each simple zero of Q. The 
three trajectories through z\ emanate from z\ at angles 9 that satisfy 

39 = n-Q\ Zl ) (mod27r). 

From the explicit formula of Q and z\ we find Q'(zi) = —y/2 — Ai and the 
three angles at Z\ are 

1 , IS 2/C7T 

= --arctan(2V2) + — -, fc = 0,l,2. 
o o 

We let 7 be the trajectory that emanates from z\ at angle 

9 = — arctan(2v / 2) = -0.4103 • • • . 
3 

D{z) = —M*) =—. f Q l/ \s)ds, zeC\ 1 . 

TXl TCI J Zl 

Figure H] shows the level curves KeD(z) = 0, KeD(z) = 1, and 1m D(z) = 
0. In order to prove that z\ and z 2 are indeed connected by 7 (as suggested 
by the figure), we use arclength parametrization of 7 

7 : z = z(t), z(0) = z\. 

Then 

z(t) 

Q l l 2 (s)ds = nif{t), 

*(o) 

where f(t) is real. Differentiating and squaring, we obtain 

Q(z(t))[z'(t)] 2 = -Af'(t)] 2 , 



Let 
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Figure 4: Contour lines of D{z) in the cubic case. In solid line \m.D(z) = 0, 
and in dashed line, ReD(z) = and ReD(z) = (left and right of the figure 
respectively). 

which implies that 

arg Q(z(t)) +2argz'(t) = n (mod27r). 

The level lines ReQ(z) = and ImQ(z) = are shown in Fig. The 
two level lines intersect of course at the zeros zi, Z2 and —i of Q. The critical 
trajectory 7 starts at z\ at an angle #0 = — | arctan(2y2), and therefore 7 
enters the shaded region of Fig. [5], which is the region where ReQ(z) < 
and ImQ(z) < 0, hence — n < argQ(z) < —n/2. As a consequence we have 
that 

-~ < axgz'(t) < 

as long as z(t) is in the shaded region. This implies that the real part of z(t) 
increases faster than the imaginary part decreases. It follows that the part 
of 7 that is in the shaded region is contained in the triangle with vertices 
Z\ = —\/2 + i, i and (1 — \/2)i. Hence 7 leaves the shaded region at a 
point on the imaginary axis above the other critical point z — — i- Then by 
the symmetry with respect to the imaginary axis we conclude that 7 indeed 
connects z\ and zi- 

This completes the proof of Theorem 12.11 
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Figure 5: Contour lines of Q(z) in the cubic case. In solid line ReQ(z) = 0, 
and in dashed line ImQ(z) = 0. 



4 Proof of Theorem 12.2 



We start by presenting a characterization of the curve T that is equivalent 
to the S'-property. 

For complex z, we define the g-f unction 

g(z) =J log{z-s)d fJl {s), (4.1) 

which is analytic when z G C \ T. We observe that Re g(z) = — £/ M (z), 
where If 1 is the logarithmic potential as defined in (11. 9j) . The equilibrium 
properties (11.101) then translate into 

Re(-g+(z)-g-(z)+V(z))=£, 2Gsupp/i, 
Re(-g + (z)-g4z) + V(z)) > £, zeT\su Wf i. 

Let us write 7 = supp /x, then from the Cauchy-Riemann equations it 
follows that the S'-property ( 11. lip is equivalent to the property that the 
imaginary part of —g + — g- + V is locally constant on 7, that is 

lm(-g + (z) - g_(z) + V(z)) = £, z G 7 (4.3) 
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with a possibly different constant i on the different components of 7. Then 
as a consequence of (14.21) and (14. 3 p we have that 

-g + ( Z )-g_( Z )+V(z)=£ + ii (4.4) 

is constant on each connected component of 7. Differentiating (I4.4p we obtain 

-g' + ( Z )-g'_( Z ) + V'(z)=0, ze 1 . (4.5) 

Next we observe that the function ^V'(z) — g'{z) is analytic for z G C\7, 
and furthermore using (14.51) : 

- g'(z)) + = -\V{z) + g\z)- = -^V'(z) - </(*))- 

for z G 7. Hence — g'(z) has a multiplicative jump of —1 on 7, and 

therefore 

Q{z):={±V\z)-g\z)) 2 (4.6) 

is analytic in the whole complex plane. The asymptotic behavior of Q(z) for 
z — > 00 follows from the fact that V(z) is a polynomial and 

g\z) = [ —dfM(s) = - + O (\), z^oo, (4.7) 



z 



since fi is a probability measure on 7. The general case of Liouville's theorem 
implies that Q is a polynomial of degree 2r — 2 if degV^ = r. 

In general this is not enough to determine the curve 7, and we need more 
information on the roots of Q(z) (or extra assumptions, such as that we are 
in the one-cut case). Since Q 1 ^ 2 {z) is analytic in C \ 7 we can deduce that 
any zero of odd multiplicity of Q is in 7. Zeros of even multiplicity can be 
anywhere and are typically not in 7. 

Let Z\ be a zero of Q of odd multiplicity. From (j4.7p we see that for z G 7, 
in the same connected component as z\, we have 



/ Q 1 l 2 {s)ds g«I 

Jzi 



This is the condition that characterizes a trajectory of the quadratic differ- 
ential Q(z)dz 2 , emanating from a zero Z\ of Q. 
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In the case V(z) = —iz 3 /3, it follows from ( 14. 6 j) and ( 14. 7p that Q should 
be taken polynomial of degree 4 so that 

QW = (44 + °(^)) 2 = -T + " + c - < 48 > 

where C needs to be determined. In order to do this, we make the assumption 
(to be justified later) that we are in the one-cut case, that is we assume that 
7 is a single curve. The endpoints of the curve are then simple zeros of Q. 
Since Q has degree four there are two more zeros, which in the one-cut case, 
should combine into a double zero. 

In our case, there is a symmetry about the imaginary axis, and therefore 
the double root should be on the imaginary axis, say z = z , and two simple 
roots are symmetric with respect to the imaginary axis, say Z\ and z 2 = — ~%\- 
This leads to 

Q(z) = --{z - z ) 2 (z - z x ){z + z-l), 
which combined with (14. 8 j) yields 

Zq = — i, and z 1 = —V2 + i, z 2 = V2 + i- 
The free constant is C = —3/4. Therefore 

Q{z) = -~{z + i) 2 {z 2 -2iz-3), (4.9) 

and we recover ( 12. ip . 

Once we have Q(z), we may obtain [i in the following way. From ( 14. 6 p 
it follows that there is an analytic branch of Q 1 ^ 2 (z) for z E C \ 7 which 
behaves as |V(z) for large z. Choose an orientation on 7. The orientation 
induces a +- side and a — side on 7, where the +- side ( — side) is on the left 
(right) as one traverses the contour according to its orientation. 

Lemma 4.1. Given the critical trajectory 7 and the polynomial Q(z), then 

— Q+ 2 (s)ds = dix (4.10) 
is a probability measure on 7. 
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Proof. The measure fi is a priori complex, however, by the construction of 7 
we have that z 

I dn(s) = — [ Q 1 l 2 {s)ds e R 
Jz x ni J Zl 

for every z e 7, so that 11 is a real measure. 
Taking z = z 2 we can compute 

/" Z2 1 f Z2 1/2 

/ d//(s) = — / Ql {s)ds 

J Z1 m J Z1 

by contour integration. Indeed, we have 



- r Q f{ s )ds = ^- [ Q l ' 2 {s)ds 



where C is a closed contour in C \ 7 that encircles 7 once in the clockwise 
direction. Moving the contour to infinity, and using the behavior of Q 1 ^ 2 at 
infinity, see (12. 3p . we find 

/i( 7 )= [ 2 d^s) = l. (4.11) 

Then if £ € [0, 1] 1— > z = z(t) is a smooth parametrization of 7 with 
z(0) = z\ and z(l) = Z2 we have that 

1 r (t) ^ 



t e [0,1] — / Q l l 2 {s)ds 

™ J Z1 



is real valued, with values for £ = and 1 for £ — 1. The derivative 
z'(t)Q 1 / 2 (z(t)) is non-zero for < £ < 1. Therefore 

*(*) 



- f gy 2 ( s )rf s 



is strictly increasing, and it follows that \i is a probability measure. □ 

Lemma 4.2. Le£ T = 71 U 7 U 72 &e de/med as in GZM . (E5J) and (EE}, 

£/ien £/ie measure /i defined by (14. lOf) zs £/ie equilibrium measure on T in the 
external field Re V. 
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Proof. From another residue calculation, similar to the one leading to ( 14.1 ip 
and based on 02. 3p . it follows that 

/ YT- S d ^ = \ V '^ - Ql/2 z eC\j. (4.12) 
Then we have 

g(z) = / log(z - s)dfi(s) 



is such that (14.51) holds, which after integration leads to (14. 3 p and to the first 
line of (Ojl . 

We extend 7 to an unbounded contour r = 7 U 71 U 72 as in section 2.2. 
The unbounded pieces 71 and 72 are such that (I2.5P and (12.61) hold. This 
leads to the second line of (14.21) . For example, if z G 72, then by (12. 6 p and 



< 2 f Z Q 1/2 (s)ds = f\v'(s) - 2g'(s))di 

J zo J zo 



= {V{z)-2g{z))-{V{z 2 )-2g{z 2 )) 
= V(z)-2g(z)-(£ + it), 

which by taking the real part indeed leads to the inequality in (14. 2 1) . 

Because of (14.21) we have that fi is the equilibrium measure on T in the 
external field Re V. □ 

Finally, because of (14. 3p we conclude that the contour T has the S- 
property, and this completes the proof of Theorem 12.21 



5 Proof of Theorem 2.4 



5.1 Riemann— Hilbert problem 

The orthogonal polynomial P n {z) characterized by (11.61) appears as the (1, 1) 
entry of the solution Y(z) of a 2 x 2 matrix- valued Riemann-Hilbert problem, 
see [12]. 

From this Riemann-Hilbert problem, the Deift-Zhou steepest descent 
method performs several explicit and invertible transformations that allow 
us to obtain asymptotic results for the entries of the matrix Y, and in par- 
ticular for P n (z), as n — > oo uniformly in different regions of C, see [10J. In 
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the present case the analysis is quite standard, except for the fact that we 
are working on a complex curve T instead of on a part of the real line. For 
this reason, we give a brief sketch of the method and refer the reader to [10J, 
[9] and [18] for the general theory involving orthogonality with respect to 
exponential weights and also for more details on a similar problem. 

We are interested in a matrix-valued function Y : C\T — > C 2x2 such that 

• Y(z) is analytic for z G C \ Y. 

1 p -nV{z) 

Y + (z) = YM)[t , ), former, 



1 



Y(z) = (l + 0(±))^ z ° n ),as^oo. 



As before, T is the contour T = 71 U 7 U 72 consisting of the critical 
trajectory 7 and its analytic extensions 71 and j 2 - See Figure |3j 

This Riemann-Hilbert problem has a unique solution if and only if the 
monic polynomial P n {z), orthogonal with respect to the weight function 
w(z), exists uniquely. If additionally P n -i(z) exists, then the solution of 
the Riemann-Hilbert problem is given by: 



Y(z) 

where 



P n {z) (CP n w)(z) 
-2m^ x Pn-\{z) -jn-i(CP n -iw)(z] 



^-^iSM da (5 - i) 



is the Cauchy transform on T, and the coefficient 7 n _i is defined as 



7n-l 



Pn-l(s)w(s)ds 

r 



-1 



5.2 First transformation 

The first transformation Y 1— > T is a normalization at 00. We use the func- 
tions 0i and 02 as in (12.51) and (12.61) which are analytic in C \ (71 U 7) and 
C \ (72 U 7) respectively and satisfy 2 (z) — (f>i(z) = ±ni for z G C \ T. We 
set 

aa-k E ^ M r- |wi ° ln \ (5. 2 ) 
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1 

7 



Figure 6: The jump matrices for the Riemann-Hilbert problem for T. 
Now, using (I4.6p . we obtain by direct integration from (12. 6 p that 



<f> 2 {z) = \V{z) - \og{z) -l + o(^\ 
for some constant of integration I. It follows that 



z — > oo, (5.3) 



e n[<h(*)-$V(z)] = z -n e -nl [ l + [\\\ : ^ , x . 



Z t 

Hence T satisfies the following Riemann-Hilbert problem: 

• T(z) is analytic for z in C \ T; 

• T has the jumps indicated in Figure [6j 

• T(z) =I + 0(\) as z-^oo. 

5.3 Second transformation 

The second transformation of the Riemann-Hilbert problem is the so-called 
opening of lenses. From the Cauchy-Riemann equations, it is possible to 
show that the sign pattern for Re 02 is as shown in Figure El Since <\>\ = 
(f>2 ± the sign pattern for Re <pi is exactly the same. 

In the second transformation we open a lens-shaped region around 7 as 
in Fig. [HI so that the lens is contained in the region where Re 02 < 0: 
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Re 02 < 



Figure 7: The sign of Re 0i = Re 02 in various parts of the complex plane. 
The solid curves are where Re 02 = 0. The curves 71 and 72 are shown with 
dashed lines. We have that 0i is real and positive on ji and 02 is real and 
positive on 72. 



We define 



S= < 




in the upper part of the lens, 
in the lower part of the lens, 



(5.4) 



Then S satisfies the following Riemann-Hilbert problem: 

• S(z) is analytic for z G C \ Fg, where T$ consists of T plus the lips of 
the lens; 

• S has the jumps indicated in Figure EJ 

• S(z) = I + O (~) as z -> 00. 



5.4 Construction of parametrices 
5.4.1 Global parametrix 

Now we seek an approximation to S that is valid for large n. The approxima- 
tion will consist of two parts, a global parametrix N away from the endpoints 
z\ and Z2 and local parametrices P at Z\ and Z2- 

The global parametrix iV satisfies a Riemann-Hilbert problem with the 
same constant jump on 7. Then R = SN' 1 will be analytic across 7. We 
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Figure 8: The contour T$ and the jump matrices on T$ in the Riemann- 
Hilbert problem for S. 

define N as 

where is given by (12.101) , see [9] , [10] and [16] . 

Then N satisfies the following Riemann-Hilbert problem: 

• N(z) is analytic for z G C \ 7; 

• N + = iV_ J) on 7 ; 

• iV(z) = I + 0(±) as 2 -^00; 

• JV(*) = 0(\z - %|~ 1/4 ) as z -> z,- for j = 1, 2. 

It is clear that iV cannot be a good approximation to 5 near the endpoints 
of 7, since it blows up at z = Z\ and z = Z2, while S remains bounded there. 
For this reason we need a different local approximation near the endpoints. 

5.5 Local parametrix 

The local parametrix P is constructed in neighbourhoods of the endpoint 
z = Zj, j = 1,2, say 

Us(zj) = {zeC \ \z- Zj \ < 5}, J = 1,2, 
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Figure 9: The jump matrices in the Riemann-Hilbert problem for P defined 
in the neighbourhood Ug{z2) of z%. 



with some small but fixed 5 > 0. We describe here the construction of P in 
Us{z2), the construction in Ug(zi) being similar. 

The local parametrix P should satisfy the following Riemann-Hilbert 
problem: 

• P(z) is analytic for z G Us(z 2 ) \ T s with a continuous extension to 



• P has the jumps on Ts PI U$(z2) as shown in Fig. (these are the same 
jump matrices as in the RH problem for S); 

• P(z) — (/ + (^)) N(z) as n — > oo, uniformly for z G dUs{z2)] 

• P{z) remains bounded as z — > Z2. 

The construction of P is given in terms of the Airy function Ai and its 
derivative. We put 



U s (z 2 )\r s ; 



P(z) = E n (z)A(n 2 / 3 f(z)) 







e 




) 



(5.6) 



where A((), f(z), and E n (z) are described below. 
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Figure 10: Contours and jump matrices in the Riemann-Hilbert problem for 
Airy functions. 

Airy parametrix A(() The matrix-valued function A(() is the solution 
of the Airy Riemann-Hilbert problem, which is posed on four infinite rays in 
an auxiliary £-plane as follows: 

• A(C) is analytic for ( € C, argC £ {0,2vr/3, -27r/3,vr}; 

• A has the jumps on the four rays as shown in Fig. [101 

• As £ — > oo, we have 

cv^o o(-°w)(T 

(5.7) 

• remains bounded as ( — > 0. 

The solution of this Riemann-Hilbert problem is given by the Airy func- 
tion Ai(C) and rotated versions of it, see [TO] and [U Sec. 10.4]. Let 

Ito(C) = Ai(C), yi(C)=wAi(wC), 2/ 2 (C) = ^ 2 Ai(a; 2 C), 
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2ni 

where u = e 3 . These are three solutions of the Airy differential equation 
y" = (y satisfying the connection formula yo(() + yi(C) + 2/2(0 = 0- F° r 
instance, in the sector < arg£ < 27r/3 we have 

A(C) = ( V?® - y $) . (5.8) 

The solution in the other sectors is obtained from this by applying the ap- 
propriate jump matrices. 

Conformal map f(z) The map f(z) is defined by 



/(*) 



2/3 



(5.9) 



which is a conformal map in a neighourhood of z = z 2 . It is assumed that 
5 > is sufficiently small so that / is indeed a conformal map on Us(z 2 ), and 
also that the lens around 7 is opened in such a way that the lips of the lens 
inside U$(z 2 ) are mapped by £ = f(z) to the rays arg£ = ±2n/3. This can 
be done without any loss of generality. 

Analytic prefactor E n (z) The prefactor E n (z) in f)5.6p is defined by 
1/1 -i\ (n l l G f(zfl A 



E n (z) - N(zj-^= x j y n - 1/6f{z) -l/4 

1/1 -i\(nV«f{zy/^-\z) \ 

~ V2 \-i 1 J V n-ysf{z)-^f){z)) ' t5 ' iUj 

which is analytic in Us(z 2 )- It is chosen so that the matching condition 
P(z) = (I + 0(l/n))N(z) for z e dU 5 (z 2 ) is satisfied. 

Then with these definitions it can be shown that P defined by (15. 6p indeed 
satisfies the Riemann-Hilbert problem for P. 



5.6 Third transformation 

In the third and final transformation we use the global parametrix N(z) and 
the local parametrices P(z) to define 



S(z)N(z)-\ zeC\(T s UU 5 (z 1 )UUs(z 2 )), 

R(z) = < (5.11) 

S(z)P(z)-\ zE(U s (z 1 )UU s (z 2 ))\T s . 
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N 



( 




N- 1 



Figure 11: The contour T R and the jump matrices on Tr in the Riemann- 
Hilbert problem for R. 

Then R has an analytic continuation across 7 and across the parts of 
S5 that are inside the disks Ug(zi) and Us(z 2 ). It satisfies the following 
Riemann-Hilbert problem: 

• R(z) is analytic for z G C \ T R , where T R is the contour shown in 



• R has jumps on each part of T R with jump matrices as indicated in 
Fig. EU 

• R(z) = I + O (i) as z ->■ 00. 

The jump matrices in the Riemann-Hilbert problem for i? tend to the 
identity matrix as n — > 00. Indeed, since P(z) = (I + 0(l/n))N(z) as 
n — > 00, uniformly for z G dU${z\) U dUs(z 2 ), we have that 



as n — > 00. On the remaining parts of Tr we even have for some positive 
constant c > 0, 



Fig.HB 





R 



{z) = R_{z) (l + 0(e- cn )) 



zeT R \(dU s ( Zl )UdU s {z 2 )) 



26 



as n — > oo. Thus the jumps on R tend to the identity matrix uniformly, and 
in fact also in L 2 (T^). 

It then follows from the general theory, see [9] and |16j . that the solution 
to the Riemann-Hilbert problem for R exists for all large enough n with 

R(z)=I + o(^J, as n ^ oo, (5.12) 

uniformly for z G C \ Tr. 



5.7 Proof of Theorem 12.41 

Once we arrive at this result for R, it is possible to reverse all the transfor- 
mations Y i— > T \- > S i— > R, since they are all explicit and invertible. The 
first thing that follows is that the original Riemann-Hilbert problem for Y 
has a unique solution for large enough n. Since 

P n (z) = Y 11 (z), 

this proves that the orthogonal polynomials P n indeed exist for every large 
enough n. 

The asymptotic formula (15.121) for R further yields the first term in an 
asymptotic expansion of Y as n — > oo. Following the effect of the inverse 
transformations R^-S^-T^-Yon the asymptotic formula (15.121) for R, 
we obtain the asymptotics of Y and therefore of P n in the various regions of 
the complex plane. This will give the different parts of Theorem 12.41 



5.7.1 Proof of part (a) 

Let z G C \ 7. We then may and do assume that the lens around 7 and 
the neighborhoods Ug(zi) and Us(z 2 ) are chosen so that z is in the outside 
region. 

From (15.111) we have that S(z) = R(z)N(z). Also, because we are outside 
the lens, we have S(z) = T(z) from (15.41) and T(z) in terms of Y(z) follows 
from (I5.2p . Combining all this we find 
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Then, for the (1, l)-entry the first part of the theorem follows in a straight- 
forward way, since 

g{z)= l -V{z)-Hz)-l- 
5.7.2 Proof of part (b) 

For z inside the lens, but outside of the two disks, we have S(z) = R(z)N(z) 
as before. From (15. 4p we then get 

n*)=s(*)( ±e U J) = *c*w*) ( ± jL J) 

where the + sign (— sign) is taken in the upper (lower) part of the lens. 
Using (15.21) and ( 15. 12ft we then find 

/ e -m Q \ / / 1 \ \ / e -n[Mz)-h v (4 o \ 

Then for the (1, l)-entry we obtain from this 

P n {z) = Y n (z) = e»[4M (l + O (J) , JV(s) 

-i] J e -n^Wjv n ( 2 ) ± e n ^ 2(z) iV 12 (2) + O ( - 



as n — > oo. This proves part (b) of the theorem. 
5.8 Proof of part (c) 

In the neighbourhoods Ug(zx) and Us(z 2 ) of the endpoints Z\ and z 2 we use 
the local parametrix P(z) to obtain an approximation for P n (z) in terms of 
Airy functions. Indeed, by (15. lip and (15.121) . 

S(z)=R(z)P(z) = (l + O^P(z) 

for z G Us(zi) U Us(z 2 ). If we assume that z is inside the disk Us(z 2 ) but 
outside the lens around 7, then we find by following the transformations (15.4)) 
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and fl52]) that 



Y(z) = n _ n , 1 1 + O - P(z 



\nJJ KJ \ e^M-^l 

Using (15.61) and (15.101) . we obtain from this that 



n- 1 /6/( 2 )-i/4 /3 (^ 



x A{n 2 ' z f{z)) 

To evaluate A(n 2 / ?J f(z)) we use (I5.8P and it follows that 

p n (z) = (i o)r^(j) 



V(z) 

e n ^ 



e~ n ^ 



^e^](l + (I), 0(1))^. , 



1 -z 



X 



n 1 /6/( z )i/4 /3 ( z )-i o \ f Ai(n 2 / 3 f(z)) 

o n- 1 / 6 /(z)- 1 /V(^)y V-* Ai V /3 /(*)) 



as n — > oo. This proves part (c) of the theorem in case z G Us(z 2 ) is outside 
the lens. A similar calculation leads to the same expression in case z is inside 
the lens. This completes the proof of part (c) of Theorem 12.41 

6 Concluding remarks 

We have presented a Riemann-Hilbert analysis of a family of polynomials 
orthogonal with respect to a varying exponential weight on certain curves of 
the complex plane. The problem was motivated by the fact that the zeros of 
these polynomials are complex Gaussian quadrature points for an oscillatory 
integral on an interval [a, b] C K. The zeros cluster on analytic arcs in the 
complex plane, which are given by a critical trajectory of a suitable quadratic 
differential. 
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We have focused on the case where the weight function is V(z) = —iz 3 /3, 
for which we were able to obtain explicit expressions throughout the Riemann- 
Hilbert analysis. A similar procedure (with more complicated computations) 
can be applied in principle to the more general case V(z) = —iz r /r with 
r > 5 and odd. The only difficulty is the determination of a curve with the 
5-property in this more general case. It would be interesting to know if we 
are in the one-cut case for every odd r. 

It is also worth remarking that the Riemann-Hilbert analysis can provide 
more detailed asymptotic information than the one given before, following the 
ideas exposed in [10], [20]. The importance of these results from a numerical 
point of view is currently under investigation. 
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